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ABSTRACT 


We discuss single top quark production via u'di — tb at the Fermilab Tevatron in the min- 
imal supersymmetric model with R-parity violation. We find that within the allowed range 
of coupling constants, the lepton-number violating couplings can give rise to observable effects 
when the slepton mass lies in a specific narrow range. For the baryon-number violating cou- 
plings, the contribution to the production rate can be quite large in the presently allowed range 
of the coupling constants. 


We show that the measurement of single top production at the 
upgraded Tevatron can be used to constrain a linear combination of products of the R-parity 
violating couplings. 


PACS: 14.65.Ha, 14.80.Ly 
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1. Introduction 


Single top quark production through W-gluon fusion[1] g+W — t+b, and quark-antiquark 
annihilation via virtual s-channel W(2] u* + dî — t +b, where i,j are the generation indices, 
are interesting to study at the Tevatron. In contrast to the QCD process of tt pair production, 
the single top production involves the electroweak interaction and can, therefore, be used to 
probe the electroweak theory. In particular, the quark-antiquark annihilation subprocess can 
be used to study models of new physics in connection with the top quark. 

It has been shown[3] that the signal for u’d’ — tb is potentially observable at the Tevatron 
with 2-3 fb”* integrated luminosity, although it will be overwhelmed at the LHC by the large 
background from tt production plus single tops from the W-gluon fusion[1]. Compared to 
Wg > tb, utdi > tb has the advantage that the cross section can be calculated reliably because 
the quark and antiquark structure functions at the relevant values of x are better known than 
the gluon structure functions entering in the calculation for the W-gluon fusion cross section. 
With a 30 fb”! data sample from Run 3 at the upgraded Tevatron with y/s = 2 TeV, it is 
possible to measure the ratio of single top production ( utd? — W — tb) and Drell-Yan cross 
section ( utdi > W — lv) to an accuracy of +8% [4]. Thus new physics effects that produce 
larger than 16% effect on the cross section ratio should be detectable[4]. 

In the Standard Model, the cross section of single top quark production ud? — tb has been 
calculated to one loop order[5]. Recently, this process has been used to study the effects of 
model-independent new physics involving the third-family quarks[6] and it was shown that at 
the upgraded Tevatron it can be a powerful probe of new physics. Effects in specific new physics 
models on this production process were evaluated in Refs.|7,8]. In the general two-Higgs-doublet 
model and the minimal supersymmetric model (MSSM) with R-parity conservation[7,9,10], the 
contribution to the production rate can reach the observable level only for tan 3 < 1 which is 
generally not regarded as viable. In most of the new physics models considered in the literature 
the enhancement effect of the single top quark production is smaller than 15% and some of the 
models can even suppress the production rate[8]. An exception to this general result can exist, 


as analysed in Ref.[8], in a non-commutting ETC model with an extra weak gauge boson of 


mass no less than 500 GeV. 

In this paper we focus on the R-parity violating MSSM[11] and evaluate the effect of R- 
parity violating couplings on single top quark production at the Tevatron. While this is an 
interesting problem in its own right, the recent anomalous events at HERA provide an additional 
motivation for the study of the R-parity violating supersymmetric couplings. The HERA data 
showed excess events in deep-inelastic positron-proton scattering at high-Q? and high x, which 
are in apparent conflict with the Standard Model expectations[12]. The excess events have 
been interpreted as evidence of R-parity breaking supersymmetry [13]. Hence the examination 
of effects of R-parity breaking supersymmetry in other processes are desirable. In Sec.2 we 
present the Lagrangian for R-parity violating couplings. In Sec.3 we evaluate the effects of 
R-parity violating couplings on the production rate of single top quark. In Sec.4 we present 


the numerical results and some discussions. 
2. Lagrangian of R-parity violating couplings 


In supersymmetric models, the R-parity of a field with spin S, baryon-number B and lepton- 
number L is defined to be 
R= (DATER, (1) 


R is +1 for all the SM particles and -1 for all super particles. R-parity invariance is often 
imposed on the Lagrangian in order to maintain the separate conservation of baryon-number 
and lepton-number. Imposition of R-parity conservation has some important consequences; 
super particles must be produced in pairs in collider experiments and the lightest super particle 
(LSP) must be absolutely stable. Thus the LSP provides a good candidate for cold dark matter. 

Despite the above mentioned attractive feature of R-parity conservation, the conservation is 
not dictated by any fundamental principle such as gauge invariance and there is no compelling 
theoretical motivation for it. The most general superpotential of the MSSM, consistent with 


SU(3) x SU(2) x U(1) gauge symmetry and supersymmetry, can be written as 
W = Wr + Wa, (2) 
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where Wp is the R-parity conserving part while Wa violates the R-parity. They are given by 


We = hiyL HE; + hiQ¡H2D5 + hi ¿Qi UF, (3) 


J 


Here L;(Q;) and E;(U;, D;) are the left-handed lepton (quark) doublet and lepton (quark) singlet 
chiral superfields. t, j,k are generation indices and c denotes a charge conjugate field. Hj 2 are 
the chiral superfields representing the two Higgs doublets. Note that the term u;L;H5 can be 
rotated away by a redefinitionH of the Higgs H, and that of the leptonic L; superfields[15] . 

In the R-parity violating superpotential Eq. (4), the A and \’ couplings violate lepton-number 
conservation, while the A” couplings violate baryon-number conservation. Ajj, is antisymmetric 
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in the first two indices and AV, is antisymmetric in the last two indices. While it is theoretically 


j 
possible to have both baryon-number and lepton-number violating terms in the Lagrangian, 
the non-observation of proton decay imposes very stringent conditions on their simultaneous 
presence|16]. We, therefore, assume the existence of either L-violating couplings or B-violating 
couplings, but not the coexistence of both. We calculate the effects of both types of couplings. 


In terms of the four-component Dirac notation, the Lagrangian of the A' and A” couplings 


that affect single top production at the Tevatron are given by 


Ly = Niza [Pidhdi + didkvi + (di) (Didi 
—E, dhuh, — ü dhet, — (di) (E, ul] + he, (5) 
Ly = Ny (EREL de + dd, + úh(di dh] + hic. (6) 


The terms proportional to À are not relevant to our present discussion and will not be considered 
here. As discussed in the proceeding paragraph, we assume that Ly and L~ are not present 
simultaneously. The Lagrangians given above can give rise to interesting effects in low energy 
processes. 

In the presence of R-parity violation, the phenomenology of MSSM changes considerably. 
For instance, the LSP is no longer stable and their decays can lead to interesting phenomenology 


2Such redefinition does not leave the full Lagrangian invariant[14], but it has no relevant consequences to 
our analysis. 


in collider experiments. Several authors have investigated the phenomenological implications 
of R-parity violating couplings at various colliders[17]. Constraints on the R-parity violating 
couplings have been obtained from perturbative unitarity[18,19], n — ñ oscillation[19,20], ve- 
Majorana mass|21], neutrino-less double 8 decay|22], charged current universality|23], e— u — T 
universality[23], 1, — e scattering[23], atomic parity violation[23], v, deep-inelastic scatter- 
ing|23], K-decay[24,25], r-decay[26], D-decay[26], B-decay[27-29] and Z-decay at LEP 1[30,31]. 
As was pointed out in Ref.[24], transforming the Lagrangian from the gauge basis to the quark 
mass basis can lead to a flavor changing neutral current (FCNC) in the up or the down quark 
sector even under the assumption of one R-parity violating coupling. FCNC processes can 


therefore provide stringent constraints on the R-parity violating couplings. 


3. Production rate of single top quark in R-parity violating MSSM 


The Feynman diagrams for the single top quark production at the tree level parton process, 
u (pi) + d (pa) — t(p3) + b(pa), are shown in Fig.1. The diagrams for the L-violating couplings 
and B-violating couplings are shown in Fig.1(b) and Fig.1(c), respectively. We assume that 
figures 1(b) and 1(c) do not exist simultaneously. Note that the B-violating couplings can also 
give rise to a s-channel diagram ut + dí > dk, > t +b. But this is a different process which 
does not interfere with the SM diagram ut + dî > W+ — t +b and it is relatively suppressed 
by the sea structure functions. Therefore, its contribution is expected to be much smaller at 
the Tevatron than the t-channel diagram in Fig.1(c). 

The amplitudes for Fig.1(a,b,c) denoted by M$ 


ij ¿MA and 6M}, respectively, are given by 


My = LL pay Ptos) u(ps)y"Prv(pa), (7) 
W 
OMY = Nor a) Paes) lps) Pt) (8) 
Er 
OMY = -AA ap Proto) = Do) Prato) 
x [D(p2) Pru" (pa) — alps) Pro (o). 9) 
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Here PL r = (17%)5)/2, É = (pı — pa)’, and the sum over k = 1,2,3 is implied. K;; are the KM 
matrix elements, $ is the center-of-mass energy squared for the parton-level process, and M dt, 
and Mex are the masses of squark and slepton, respectively. 
The SM parton-level cross section at tree level is 
4 2 (2 2)2 
pa a HRA wr (25 + MP), (10) 
where we have neglected the masses of the bottom quark and the initial partons E) The contri- 


bution of L-violating couplings to the parton-level cross section is given by 


2 
1 (ê — MP E kij k33 (11) 


aN 
AG, = - IA , 
L 


where Da is the width of the charged slepton é%. In the R-parity conserving MSSM, the 
charged sleptons € will decay into charginos and neutralinos via the processes č% — va + a 
(j = 1,2) and € > +38 (j = 1,2,3, 4), where XY; and X9 represent a chargino and neutralino, 
respectively [32]. However, in the R-parity violating MSSM, the slepton can also have the decay 
modes é% — w + d' (i, j = 1,2,3) via the A couplings. The partial widths are given by 


2 2 
a z+ g 2 2 472 
D(ez > Vek + Xj ) = 167M3, [Usa] (ma M2. ) > (12) 
zk k <0 g’ / 1 1 2 / 2 Be 
Me, — e+ Ro = SM3 sw; + E sw) No Ma — Mio > (13) 
ek 
~k =j i o ( hi)” 2 ga 
Des >u + d ) = 167 M3 (Mé = M2) > (14) 
eh 


where sw = sin Oy», cy = cos Oy and the masses of the lepton e* and down-type quark d’ 
are neglected. The masses of charginos and neutralinos, and the matrix elements U;; and Ni; 
depend on the SUSY parameters M2, Mı, u, and tan [10]. Here, Ma and M, are the masses 
of gauginos corresponding to SU(2) and U(1), respectively. u is the coefficient of the H¡ Ho» 
mixing term in the superpotential, and tan 8 = v2/v, is the ratio of the vacuum expectation 
values of the two Higgs doublets. 


3We neglected the contribution of third-family sea quark in the initial states 


The contribution of the B-violating couplings to the cross section is given by 


n E $ — M? 
Age = Ir. NE TA $+ M2 — 2M2 
Oi; 247 13k“3jk 7 32(3 — M2.) s+ t dk, 
5 M5, — M$ 8+ M3, — M$ 
2Mi—E log —— E |, 15 
E o Ma (1a) 


where the sum over k = 1,2,3 is implied. 


For the B-violating couplings, we keep only the interference term ôM, Mi . But for the L- 


violating couplings, we kept the higher order term |6M)/|? since the interference term 6M. Mi ~ 


Majma which drops out when the mass of d’ is neglected. 


The total hadronic cross section for the production of single top quark is obtained by 


Aan a aca E (16) 


ij 7 TS dt 


where s is center-of-mass energy squared, 7 = (M; + M,)?/s and 7 is defined by T = 2122 with 
X12 denoting the longitudinal momentum fractions of the initial partons ¿ and j, respectively. 


The quantity dL;;/dr is the parton luminosity defined by 


di [Aen 1012 (0/20) + (A B) (17 


dt Em 
where A and B denote the incident hadrons and the functions fA and f are the usual parton 


distributions. 


5. Numerical results and conclusion 


In our numerical calculation, we use the CTEQ3L parton distribution functions[33] with 
u = V8, and assuming M, = 175 GeV and vs = 2 TeV. The KM matrix elements Ky = 


— Kə = 0.22 are used in our calculation. 


5.1 L-violating couplings 


For the contribution of L-violating couplings, we assume the masses of sleptons é% to be 
degenerate. The contribution of L-violating couplings to the cross section is sensitive to the 
width of the charged slepton ef, which depends on the SUSY parameters Mə, Mı, p, tan 6, 
and all Y couplings. 

In our calculation we use the GUT relation M, = LM; N 35M». Since the masses of 
charginos and neutralinos are not sensitive to tan 5, we fix tan 5 = 2 and retain M and y as 


variables. 


The upper limits of the L-violating couplings for the squark mass of 100 GeV are given by 


ll E OA RSS LO, (18) 
Xy OG GEL; (19) 
aa ee. OE (20) 
Ae) = O16, G23), (21) 
ll Ube HA 2S), (22) 


The first set of constraints come from the decay K — mvv with FCNC processes in the down 
quark sector[24]. The second and forth set of constraints are obtained from the semileptonic 
decays of B-meson[29]. The third constraint, i.e., that on the coupling Ms, is obtained from 
the Majorana mass that the coupling can generate for the electron type neutrino|21]. The last 
set of limits have been derived from the leptonic decay modes of the Z[30]. 


Also the following constraints are derived for the combinations of A” couplings 


MM; Ag3Aga < 11x 107°, (6 = 1,2; 7 = 1,2, 3), 23 
AiimAjano Min2 Aint < imãs (i jn = 1, 2,3), 24 


AA AA» Alo Ago < 107%, 
Nios Acero. Asia Aion < 1077, 
where the first set of constraints are derived in Ref.|29] and the other derived in Ref.[25]. 


Inputing the upper limits of the relevant L-violating couplings for squark mass of 100 GeV, 


we get the maximum contribution of L-violating couplings which are shown in Fig.2 and Fig.3. 
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Figure 2 shows the histogram of the differential cross section versus the invariant mass of the tb 
system over a bin size of 10 GeV with the favorable parameters Mz = —u = 200 GeV. The solid 
line is for the standard model. To illustrate the contributions of sleptons of different masses, we 
superpose the effects of three sleptons on the same curve. The dashed, dotted and dash-dotted 
lines are the standard model plus the slepton contributions for three different slepton masses 
respectively: 230 GeV, 300 GeV and 350 GeV. The resonance behavior is already manifested. 
Because of their narrow widths, for each slepton the contributions of the A'-couplings are 
negligible for a couple of bins away from the resonance. This will help to identify the signal of 
the slepton production. 

Figure 3 shows the ratio of the total slepton contribution integrated over 3 to that of the 
standard model as a function of the slepton mass. For slepton mass lighter than 180 GeV, the 
contribution is very small. When the slepton mass becomes heavier than 180 GeV, which is the 
threshold for the slepton to decay into t and b, the contribution increases sharply and reaches 
its maximum sizes at the slepton mass of about 230 GeV. Then the effect drops quickly with 
further increase in the slepton mass. The large slepton mass suppression can be understood as 
follows: When the slepton mass is large the parton cross section contributions coming mainly 
from $ ~mass of slepton require large momenta from the initial partons which is suppressed 
by their structure functions. An additional suppression is caused by the increase of the slepton 
width when the slepton mass increase. 

Figure 3 also shows that the contribution to the production rate increases with the increase 
of Mə and u. This is because with the increase of Mə and u, the chargino and neutralino 
masses increase and thus the width of the slepton decreases. The overall contribution to the 
production rate can exceed 20% when slepton mass lies in a narrow range ( 200 GeV~ 270 
GeV) and M> = —y > 200 GeV. As shown in Fig.2, the contribution is confined within a bin 
of 10 ~ 20GeV. So, one can use the upgraded Tevatron to search for enhanced single top quark 


production so as to further constrain the L-violating couplings in this mass range of sleptons. 


5.2 B-violating couplings 


In our calculation, we assume the masses of dt, to be degenerate. Then the contribution of 


B-violating couplings can be parametrized as 


Aoyn 
E E + Fan Azas Aga A (27) 
where F;; depend on squark mass and are shown in Fig.4. From Fig.4 we see that Fy, and Fy, 
are the same and negligibly small. Since Ff is smaller than FY, and A/,3 < 107*[19,20], we can 
also neglect the term F141134310. 
For squark mass of 100 GeV, we take the upper limits for other B-violating couplings as 


given by[18,19,31] and are derived from perturbative unitarity and Z decay. 


312 < 0.97, (28) 


i325 3035 A213 < 1.25, (29) 


From these upper limits, we find the upper limit of the contribution to be Ao" /o? < 2752%. 
This shows that the contribution can be quite large in the allowed region of B-violating cou- 
plings. Of course, in this upper limit of couplings, our results are no longer reliable and the 
higher order terms such as (A(3244,2)? have to be included. 

Neglecting the terms proportional to Hj, Fs, and F; we present in Fig.5 the plot corre- 
sponding to Ao /c9 = 20 % in the (A324317, Mg) plane. The region above the plot corresponds 
to Ao /09 > 20 % while the region below the plot corresponds to Ao /a9 < 20 %. From Fig.5 we 
see that if we assume an observable level of 20% at the upgraded Tevatron, the coupling A{s2A519 
can be probed down to 0.01, 0.02, 0.03, 0.04, 0.06, 0.08, 0.1 and 0.13 for Ma =100, 200, 300, 
400, 500, 600, 700 and 800 GeV, respectively. As mentioned above, the present upper limit is 

1324319 < 1.25 x 0.97 ~ 1.2 for squark mass of 100 GeV. So, the single top production process 
at the upgraded Tevatron can be used meaningfully to probe the product of the B-violating 


couplings \",, and Af... 
132 312 


In summary, we studied the single top quark production via ud — tb at the Fermilab 
Tevatron in R-parity violating supersymmetry. We found that within the allowed range of the 


coupling constants, the effects of the L-violating A’ couplings can be observed at the upgraded 


Tevatron only for slepton mass in a narrow range. As shown in Fig.2 a distinctive resonance 
is associated with the enhanced production of single top. This can also serve as a signal 
for the production of slepton via L-violating couplings or further constraint to their coupling 
strengths. For B-violating A” couplings, in the allowed range of the relevant coupling constants 
Eqs. (23/29), the contribution to the production rate can reach the observable level for a wide 
range in squark mass. So the upgraded Tevatron can make a powerful probe for the product of 
the B-violating couplings As, and A3,.. Failure to observe a signal of enhancement to the single 
top production will be an indication of small B-violating couplings as shown in Fig.5. We note 
that owing to the linear vs quadratic dependence of the products of R-parity violating couplings, 
the B-violating couplings will be a more sensitive probe than the L-violating couplings. 

It should be noted that the sparticle exchange effect in single top quark production ex- 
plored in the present article is complementary to the approach of direct sparticles production 
at Tevatron. This latter approach can also yield a mass limit on the sparticle. In the case of the 
L-violating coupling, a lower limit of 100 GeV on the mass of squark/gluino has been obtained 
in [34] and the limit can possibly be pushed up to 250 GeV for the accumulated dilepton data 
from RUN 1 [35]. But this sparticle production with the L-violating coupling does not allow a 
direct measurement of the slepton mass. If we use the MSSM with grand unification, in which 
the squark and slepton are related, a similar limit for the slepton mass can also be obtained 
as for the squark, i.e., a 100/250 GeV mass bound can be set. The slepton pair production 
in R-parity conserving MSSM is also a limited probe of the slepton mass at the Tevatron [36]. 
However, our approach allows a direct reach of the slepton mass, albeit a limited range of val- 
ues. For the B-violating SUSY model corresponding to the A”, the Tevatron reach of squark 
mass by direct production of the squark is very limited. The approach presented in this article 
can probe a wide range of the the ratio of the product of the A” and the squark mass. 

We conclude by noting that the s-channel squark intermediate process, u’+d? > dk, > t+b, 
is small at the Tevatron as pointed out early. But this process may not be small in the 
environment of pp collision of LHC. This process has an interesting distinctive signal of double 


b production and is currently under investigation. 
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Figure Captions 


Fig. 1 Feynman diagrams for single top quark production without associated jet other than 
the b: (a) tree-level in the SM, (b) contribution of L-violating couplings, (c) contribution of 
B-violating couplings. The blobs denote the R-parity violating SUSY vertices. 

Fig. 2 The histogram of the differential cross section versus the invariant mass of the tb 
system over the bin size of 10 GeV with the parameters Mz = —u = 200 GeV. The solid line is 
for the standard model. The dashed, dotted and dash-dotted lines are the standard model plus 
the slepton contributions for three different slepton masses respectively: 230 GeV, 300 GeV 
and 350 GeV. The vertical scale is in pb/GeV. 

Fig. 3 The integrated contribution of L-violating couplings to the cross section versus the 
slepton mass. 

Fig. 4 The form factors F; as functions of squark mass. 

Fig. 5 The value of A(3243,2 versus squark mass for Ao/oy = 20 %. The region above the 


plot corresponds to Ao /09 > 20 % while the region below to Ag/ao < 20 %. 
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